A- symmetries of Dynamical Systems, 
Hamiltonian and Lagrangian equations 

Giampaolo Cicogna 
Dipartimento di Fisica, Universitd di Pisa 

and INFN, Sezione di Pisa, 
Largo B. Pontecorvo 3, 50127 Pisa (Italy) 
cicognaOdf .unipi . it 

Abstract 

After a brief survey of the definition and the properties of A-symmetries 
in the general context of dynamical systems, the notion of "A-constant 
of motion" for Hamiltonian equations is introduced. If the Hamiltonian 
problem is derived from a A-invariant Lagrangian, it is shown how the 
Lagrangian A-invariance can be transferred into the Hamiltonian context 
and shown that the Hamiltonian equations turn out to be A-symmetric. 
Finally, the "partial" (Lagrangian) reduction of the Euler-Lagrange equa- 
tions is compared with the reduction obtained for the corresponding Hamil- 
tonian equations. 

1 Introduction (A-symmetries) 

Let me briefly recall for the reader's convenience the basic definition of A- 
symmetry (with lower case A), originally introduced by C. Muriel and J.L. 
Romero in 2001 [3 [5]. 

Let me consider the simplest case of a single ODE A(t, u{t),u, ii ...)= for 
the unknown function u = u(t) (I will denote by t the independent variable, 
with the only exception of the final Section 4, because the applications 1 am 
going to propose will concern the case of Dynamical Systems (DS), where the 
independent variable is precisely the time t, and ii = du/dt, etc.). Given a 
vector field 

the idea is to suitably modify its prolongation rules. The first A-prolongation 
aI^-* is the defined by 

aW = aW + A(^-™)^ (1) 
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where A = \{u,u,t) is a C°° function, and X'^^^ is the standard first prolonga- 
tion. Other modifications have to be introduced for higher prolongations, but 
in the present paper I need only just the first one. 

An n-th order ODE A = is said to be A-invariant under X^"'' if 



= 

A=0 



where is the n-th A-prolongation of X. 

It should be emphasized that A-symmetries are not properly symmetries, 
because they do not transform in general solutions of a A-invariant equation 
into solutions, nevertheless they share with standard Lie point-symmetries some 
important properties, namely: if an equation is A-invariant, then 

• the order of the equation can be lowered by one 

• invariant solutions can be found (notice that conditional symmetries do the 
same, but A-symmetries are clearly not conditional symmetries) 

• convenient new ("symmetry adapted") variables can be suggested. 

In the context of DS, which is the main object of this paper, the first two 
properties are not effective, the third one is instead one of my starting points. 

Before considering the role of A-symmetries in DS, let me recall that many 
applications and extensions of this notion have been proposed in these 10 years: 
these include extensions to systems of ODE's, to PDE's, applications to varia- 
tional principles and Noether-type theorems, the analysis of their connections 
with nonlocal symmetries, with symmetries of exponential type, with hidden, or 
"lost" symmetries, with potential, telescopic symmetries as well. Other inves- 
tigations concern their deep geometrical interpretation, with the introduction 
of a suitable notion of deformed Lie differential operators, the study of their 
dynamical effects in terms of changes of reference frames, and so on. Only the 
papers more directly involved with the argument considered in this paper will 
be quoted; for a fairly complete list of references see e.g. [3l|4l[5]. A very recent 
application concerns discrete difference equations [6]. 



2 A-symmetries for DS 

I am going to consider the case of dynamical systems, i.e. systems of first-order 
ODE's 

Ua = fa{u,t) (a = 1, . . . ,m) 

for the m > 1 unknowns Ua ~ Ua{t). 

Let me start with a trivial (but significant) case: if the DS admits a rotation 
symmetry, then it is completely natural to introduce as new variables the radius 
r and the angle 9, and the DS immediately takes a simplified form, as well known. 
However, in general, symmetries of DS may be very singular, and/or difficult to 
detect. An example can be useful: the DS 

ill — UlU2 U2 = —u\ 
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admits the (not very useful or illuminating) symmetry generated by (with = 
ul + ul) 

^ = [— 3- — — + 3 log — — 7r~ ■ 

V r'^ U2 + r / oui \ r-^ U2 + r/ OU2 

In this example the rotation (with a commonly accepted abuse of language, the 
same symbol X denotes both the symmetry and its Lie generator) 

d d 

X = U21^ Ul- 

aui OU2 

is a \- symmetry (its precise definition will be given in a moment), and not a 
symmetry in the "standard" sense; nevertheless^ still introducing the variables 
as before, i.e. r and 0, the DS takes the very simple form 

r — Q 9 = — rcos 6 . 

This is just a first, simple example of the possible role of A-symmetries in the 
context of DS. 



2.1 A-symmetries of general DS 

The natural way to extend the definition ([T]) of the first A-prolongation of the 
vector field 

d d 
X = ipa{u,t)-^ + T{u,t)— = ip-Vu + rdt 

to the case of m > 1 variables Ua is the following (sum over repeated indices) 

Xj^^ = X<-'^ +Aab{ipb-TUh)-yu^ 

where now A = A{t,Ua,Ua) is a (m x m) matrix; accordingly, I denote by the 
upper case A these symmetries in this context. 

To simplify, let me assume from now on t = (or use evolutionary vector 
field, it is not restrictive). 

Then the given DS is A-invariant under X (or AT is a A-symmetry for the 
DS), i.e. X'^\u - f)\u=f = 0, if and only if 

[f,f]a+dtiPa^-{ALp)a (a=l,...,m) 

where 

Given X , we now introduce the following new m+1 "canonical" (or symmetry- 
adapted) variables {notice that they are independent of A): precisely, m — 1 
variables Wj = Wj{u) which, together with the time t, are X- invariant: 

XWj^Xt^O (j 1, . . . ,TO - I) 
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and the coordinate z, "rectifying" the action of X, i.e. 

^-^^ 

oz 

Writing the given DS in these new variables, we obtain a "reduced" form of the 
DS, as stated by the following theorem [TIIHIH]. 

Theorem 1 Let X be a K- symmetry for a given DS; once the DS is written in 
terms of the new variables u)j,z,t, i.e. 

Wj = Wj{w, z,t) z — Z{w,z,t) 

the dependence on z of the r.h.s. Wj , Z is controlled by the formulas 

dWj dwj , , , , , dZ dz , , . 

One has: 

• // A = then Mj = Mm = and Wj , Z are independent of z 

• I f A = XI then only Z depends on z 

• Otherwise, a "partial" reduction is obtained: If some Mk = 0, then Wk is 
independent of z. In terms of the new variables, the K-prolongation becomes 

(^\ d d d 

oz owj oz 

The first case (A = 0) clearly means that X is an exact, or standard Lie 
point-symmetry [7] ; the second one has been considered in detail by Muriel and 
Romero [5] (notice that actually it would be enough to require Aip = Xip); the 
last case has been dealt with in [5]: several situations can be met, depending on 
the number of vanishing Mj (e.g., one may obtain triangular DS, or similar). 



2.2 Hamiltonian DS 

I now consider the special case in which the DS is a Hamiltonian DS. Obvious 
changes in the notations can be introduced: the m variables u = Ua(t) are 
replaced by the m = 2n variables qa{t),Pa{t) (a = 1, . . . , n), and the DS is now 
the system of the Hamiltonian equations of motions for the given Hamiltonian 
H^H{q,p,t): 

u = JVH = F{u,t) , V EE V„ = (V,,Vp) 
where J is the standard symplectic matrix 
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A vector field X can be written accordingly (with a = 1, . . . , 2n ; a = 1, . . . , n) 

d d 

X = (pa{u,t)- \-llJa{u,t)— = ^-\/u , ^ = i(pa,tpa) 

and all the above discussion clearly holds if X is a A-symmctry for an Hamil- 
tonian DS. Clearly, here A is a (2n x 2n) matrix. But Hamiltonian problems 
possess certainly a richer structure with respect to general DS, which deserves 
to be exploited; a first instance is clearly provided by the notion of conservation 
rules, with its related topics. 

Let me then distinguish two cases: 

(i) X admits a generating function G{u, t) (then X is often called a "Hamiltonian 
symmetry"): 

* = JVG i.e. ip = VpG , V = -VgG (2) 

this implies VD^G = 0, where is the total derivative, i.e. G is a constant of 
motion, D^G = 0, possibly apart from an additional time-dependent term, as 
well known. 

(ii) X does not admit a generating function: also in this case, defining 

S{u,t) = V • * one has that DtS = (3) 

and therefore, if S* 7^ const, then S is a first integral (the examples known to me 
of first integrals of this form are rather tricky, being usually obtained multiplying 
symmetries by first integrals; but they "in principle" exist, and their presence 
will be important for the following discussion, see subsect. 3.4). 

Direct calculations can show the following: 

Theorem 2 // the Hamiltonian equations of motion admit a A-symmetry X 
with a matrix A, then: 
in case (i) 

V(AG) = JA* = JAJVG 

in case {ii) 

DtS = -V(A*) . 
When this happens, G (resp. S) will be called a ^^A-constant of motion" . 

If A = 0, i.e. when X is a "standard" (or "exact" ) symmetry, the above 
equations become clearly the usual conservation rules; A-symmetries can then 
be viewed as "perturbations" of the exact symmetries. More explicitly, the 
equations in Theorem 2 state the precise "deviation" from the conservation of 
G (resp. of S) due to the fact that the invariance under X is "broken" by the 
presence of a nonzero matrix A. 

As a special case for case (i), the following Corollary may be of interest: 
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Corollary 1 Under mild assumptions (A^ = X^, A = X{G)), then the A- 
constant of motion G satisfies a "completely separated equation" , involving only 
G{t): 

G = i{t,G) . 

This equation expresses how much the conservation of G{t) is "violated" along 
the time evolution. If A is in some sense "small" , then G is "almost" conserved. 



3 When a A-symmetry of the Hamiltonian equa- 
tions is inherited by a A-invariant Lagrangian 

3.1 A-invariant Lagrangians, Noether theorem and 
A-conservation rules 

Let me consider (for simplicity) only first-order Lagrangians: 

C = C{qaAa,t) = 
Such a Lagrangian is said to be A'^'^^-invariant pTH [TT| under 

if there is an [n x n) matrix 

^(C) ^ A(^)((7,g,i) 

such that 

where {^X^f'^^ is the first A^'^-'-prolongation of X'^'~^ (the notation is rather 
heavy, to carefully distinguish the Lagrangian case from the Hamiltonian one, 
to be considered in the next subsection). 
We then have [TT] 

Theorem 3 // the Lagrangian C is K'''^^ -invariant under X^'^'^ then, putting 

VaP = VaPp With Pf3 = -— 



has 



AP = -Ai'=>,^ = -(A(-)^)^p. 



where P — Tr('P) = ^PaPa,' or also, introducing a "deformed derivative" Dt 
{Dt)o.p = DtSc^p + aI^^ then Tr(AP) = . 
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This result can be called the "Noether K^''^ -conservation rule". Indeed, if A'^'^^ = 
0, the standard Noether theorem is recovered. 

In the special case A^'^'^s — X(p, the above result becomes 

DtP = where Dt ^ Dt + X . 

Theorem 3 can be extended [llj to divergence symmetries and to generalized 
symmetries as well. Also, higher-order Lagrangians can be included: the A^'^^- 
conservation rule has the same form, but Vap is different: for instance, for 
second-order Lagrangians one has 

djC . . dC dC 

3.2 Prom Lagrangians to Hamiltonians 

Assume to have a Lagrangian which is A^'^'-invariant under a vector field 

d 



and introduce the corresponding Hamiltonian H with its Hamiltonian equations 
of motion. The natural question is whether the A'^-^^-symmetry X'^'-''> of the 
Lagrangian is transferred to some A^^^-symmetry AT^^^ of the Hamiltonian 
equations of motion. Two problems then arise: i) to extend the vector field 
X'^'--'^ to a suitable vector field X''-^\ and ii) to extend the (n x n) matrix A^-^^ 
to a suitable (2n x 2n) matrix A^^^^. 

First, the vector field X^^' is expected to have the form 

X^X(^) - (4) 

Oqa OPa 

where the coefficient functions ip must be determined. This can be done observ- 
ing that the variables p are related to q (and then the first A'^-^^-prolongation 
of X'^^'> is needed, where the "effect" of A*^-^) is present). One finds, after some 
explicit calculations, 

d / ^ ,(C) djC\ dA^f^ dC ding 

But the term in parenthesis vanishes if the Lagrangian is A*^ ■''•'-invariant, thanks 
to Theorem 3; in addition, if A'^'' does not depend on q (as happens in most 
cases, otherwise a separate treatment is needed, see subsect. 3.4), then we are 
left with 

^ "^^ ■ 

This implies that X admits a generating function, which is just 

G = faPa = P 



7 



using the notations introduced in Theorem 3. 

Second, let me now introduce the following (2n x 2n) matrix 

where A^^^ must satisfy (A is not uniquely defined, as well known) 

(2) dipj_ _ ,(c)d(pi 

It is well known that Euler-Lagrange equations coming from a A^^^-invariant 

Lagrangian do not exhibit in general K- symmetry. In contrast with this, it is 
not difficult to verify explicitly that the Hamiltonian equations of motion turn 
out to be A^^^-symmetric under the vector field X'^^^ obtained according to the 
above prescriptions. 

In conclusion, I have shown the following 

Theorem 4 If C is a A-invariant Lagrangian under a vector field X^^^ with 
a matrix A^^^ (not depending on q), one can extend X^^^ to a vector field 
X = X^^^ and the (n x n) matrix A^^^ to a {2n x 2n) matrix A = A^^^ in 
such a way that the resulting Hamiltonian equations of motion are Asymmetric 
under X; in addition, G = ipaPa is a A- constant of motion. 



Example 1 The Lagrangian (with n = 2) 
2 \qi 

is A^^^ -invariant under 



2 I 

+ 2 - 91 «2)^ exp(-252) + qi exp(-g2) 



d 



d 



with 



dqi dq2 
A('^) = diag . 



It is easy to write the Hamiltonian equations of motion and to check that they 
are indeed Asymmetric under 



X 



a d 

oqi dq2 



■Pi 



d 

dpi 



with 



A = A(^) 



0\ 

91 
-P2 qi 
0/ 

X-invariant coordinates are wi ~ qi cxp(— (72), W2 = qiPi, W3, ~ p2, and, as 
expected, the generating function G = W2 + W3 satisfies the A- conservation rule 



( 



-V\ 
\ 



V„AG = JA# or AG 



-9iG . 
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A special, but rather common, case is described by the following; 



Corollary 2 // 

A^^V = c(p 

where c is a constant, then also A$ = c$ and the "most complete" reduction 
of the Hamiltonian equations of motion is obtained: 

G = "fiG,t) Wj=Wj{w,G,t) z^Z{w,G,z,t) 

3.3 Reduction of the Euler-Lagrange equations versus 
the Hamiltonian equations 

In this section I want to compare the reduction procedure which is provided 
by the presence of a A-symmetry of a Lagrangian (i.e. the reduction of Euler- 
Lagrange equations) with the analogous reduction of the Hamiltonian equations 
of motion. 

Let me start recalling that any vector field X = Lpad/dqa admits n (0-order) 
invariants (as already said, see subsect. 2.1) 

''^3 ~ "^jili (j = Ij • • • J ^1 ^ 1) Ei-nd the time t 

and n first-order differential invariants rja — rja{q, t, q) under the first prolonga- 
tion 

X'^^^rja = (a = l,...,n) . 

Both if X*^^^ is standard and if it is a A prolongation {under the condition 
A (/? = A (/j) , it is well known that Wj are n — 1 first-order differential invariants 
(notice that this is an "algebraic" property, not related to dynamics). If one 
now chooses another independent first-order differential invariant C — (^{q,t,q), 
then one has that any first-order A'^'-invariant Lagrangian is a function of the 
above 2n invariants 

t, Wj,Wj and Q . 

Writing the Lagrangian in terms of these variables, the Euler-Lagrange equation 
for C, is then simply 



This first-order equation provides in general a "partial" reduction, i.e., it pro- 
duces only particular solutions, even considering the Euler-Lagrange equations 
for the other variables [TOl [3] (notice that this is true both for exactly invariant 
and for A'^'^^-invariant Lagrangians). 

I want to emphasize that, introducing A-symmetric Hamiltonian equations 
of motion along the lines stated in Theorem 4, then a "better" reduction is 
obtained, and no solution is lost. The following example clarifies this point. 
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Example 2 The Lagrangian (n = 2) 

2\qi ) 2\qi q-^J 

is K'^'''' -invariant under 



oqi dq2 



with A(^) = diag (1,1). With 

w = qiq2, w = qiq2 + qiq2, C = loggi 

the Lagrangian becomes 

and the Euler- Lagrange equation for ( is 

dC/dC, = C,=Q or g'i=giloggi 
with the particular solution 

q\{t) = exp(ce*) . 

The corresponding Hamiltonian equations of motion are Asymmetric under 

d d d d 

X = qi- q2^ pi^ ^P2^ 

oqi oq2 dpi op2 

with A = diag (1, 1, 1, 1). Invariants under this X are 

wi = qiq2, W2 = qiPi, W3 = 92^2 

and X is generated by G = W2 — W3. A "complete" reduction is obtained: with 
z = loggi, we get 

Wi = WiWs W2 = Ws — W2 

G = —G z = z + W2 — W3 
The above "partial" (Lagrangian) solution C = corresponds to 

z = z, W2 = W3 = c = const, wi = cwi . 

From the Hamiltonian equations, instead, e.g.: 

qi{t) = exp(ce*) + ci exp(— etc. 

The reader can easily complete the calculations. 
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3.4 When A^^^ depends on q 



If A depends also on q (see eq.s (|4l5p ). the calculations performed in subsect. 
3.2 cannot be repeated, the coefficient functions ipa cannot be expressed in the 
simple form ^ and the vector field X does not admit a generating function 
G. In this case one can resort to the other quantity S, introduced in (jSj, which 
provides a A-constant of motion. An example can completely illustrate this 
situation. 

Example 3 (n = 1) 

IS -invariant under 

X'-^'^ = q-^ with A(^) = q + q . 
oq 

One finds tp = ~qp — p and the resulting vector field 

does not admit a generating function. Nevertheless, the Hamiltonian equations 
of motion are A-symmetric under X with 



A 

Here 



q + q 

-p q + q 



s = -q 

satisfies DtS = — V(A$) and is a A-constant of motion. 



4 A digression: general A-invariant Lagrangians 

The A-invariance of a Lagrangian C = C{q,q,t) considered in subsect. 3.1 is 
a special case of a much more general situation. Instead of n time-dependent 
quantities qa{t), let me consider now n "fields" 

Ua{xi) (a = l,...,n; i = l,...,s) 

depending on s > 1 real variables Xi. Now, the Euler-Lagrange equations be- 
come a system of PDE's, and the notion of /i-symmetry [12l US] extends and 
replaces that of A-symmetry (or A-symmetry if n > 1). 

In this case, there are s > 1 matrices A^ {n x n), which must satisfy the 
compatibility condition 

AAj - DjK, + [A„ Kj] = (A = D.^) (7) 
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which can be rewritten putting Di = Di 5 + Ki (or, in explicit form: {Di)^^ = 
DiSa/3 + {Ai)a[3, with a notation extending the one introduced in Theorem 3), 

[D,,D,] = . 

Then one has [H [H [13] : 

Theorem 5 Given s > 1 matrices Aj satisfying Qj, there exists (locally) a 
(n X n) nonsingular matrix T such that 

A, = r-\D,r) . 

If a Lagrangian C is K-invariant under a vector field 

d 

X = ipa^ 

OUa 

then there is a matrix-valued vector 

which is A-conserved; this A- conservation law holds in the form 

Tr [r-iA(rp,;)] = 

or in the equivalent forms 

DiP^ = -(AOa/3(P.O/3a = -Tr{A^V^) , where P^ = {V,)c.c. = TrP, , 

Tr(5, r^) = . 

For first-order Lagrangians the A-conserved "current density vector" Vi is given 

by 

[Pijap = '■Pa^ where up^i = — — 

and for second-order Lagrangians by 

dC ^ 9£ dC 



Example 4 Letn = s = 2. Writing for ease of notation, x,y instead ofxi, X2, 
and u = u{x,y), v = v{x,y) instead o/mi, U2, consider the vector field 

d d 



and the two matrices 



J \Uy 
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and then 

\u 1 

It is easy to check that the Lagrangian 

^ ^ 2 ("^ ^ ^ u (""^"^ ^ "y""!/) + ""^ exp(-2w) 

is A-invariant (or better, in this context, ^-invariant) hut not invariant under 
the above vector field X. The ^-conservation law Tr(DiVi) ~ takes here the 
form 



= Dx{uux - Vx -) + Dy[uuy -Vy -) = ul+Uy 



In agreement with Theorem 5, the r.h.s. of this expression is precisely equal to 
-Tr(A,7',) = -(A,(p)„^ . 

dUaS 

Notice that in this case the quantity + is just the "symmetry-breaking 
term", i.e. the term which prevents the above Lagrangian from being exactly 
symmetric under the vector field (0). 

It should be remarked that /^-symmetries are actually strictly related to 
standard symmetries, or - more precisely - are locally gauge-equivalent to them 
(see for details [HI El III])- 

Given indeed the vector field X = ipad/dua and the s matrices A^, let me 
denote by 



J dUa,J 

the infinite A-prolongation of X, where the sum is over all multi- indices J as 
usual, and = ^pa. Introducing now the other vector field X 

~ ~ d 

X^(f>a^ with ifa = {T(p)a 

OUa 

where T is assigned in Theorem 5, and denoting by 



J dUa,J 



the standard prolongation of X, one has [131 HI] that the coefficient functions 
"^iif^ of the A prolongation of X are connected to the coefficient functions (^L'^'' 
of the standard prolongation of X by the relation 
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In the particularly simple case n = 1 (i.e., a single "field" u{xi)), then the 
s > 1 matrices Aj, and the matrix T as well, become (scalar) functions Aj and 
7; in this case, if a Lagrangian is /x-invariant under the vector field X, then 
it is also invariant under the standard symmetry X = 7X. In addition, the fj, 
conservation law can be also expressed as a standard conservation rule 

AP' = 

where = j(padC/dUa,i is the "current density vector" determined by the 
vector field X = jX. 

Example 5 Let now n = 1, s = 2, and let m,e introduce for convenience as 
independent variables the polar coordinates r,9. I am considering a single "field" 
u = u{r,9) and the rotation vector field X = d/dO. The Lagrangian 

C = ^r^ exp(-e6')w^ + ^ exp(e6')wg 

is clearly not invariant under rotation symmetry (if e ^ Q), but is ^-invariant 
with Ai = 0, A2 = e. The above Lagrangian is the Lagrangian of a perturbed 
Laplace equation, indeed the Euler- Lagrange equation is the PDE 

r'^Urr + 2rUr + exp{2e9){u0g + €u$) = . 

It is easy to check that the current density vector 

P = { — exp{—e6)urU0 , ^r^ exp(— e^)u^ 

satisfies the j2- conservation law 

DiPi = -eP2 . 

According to the above remark on the (local) equivalence of the ji-symmetry 
X to the standard symmetry X = ^X = exp{e6)d/d9, also the (standard) 
conservation law = holds, with 

^ = r^UrUe , ^r^ul - ^ eKp{2e9)ufj . 

5 Conclusions 

I have shown that the notion of A-symmetry, and the related procedures for 
studying differential equations, can be conveniently extended to the case of 
dynamical systems. 

The use and the interpretation of this notion becomes particularly relevant 
when the DS is a Hamiltonian system, and even more if the symmetry is inher- 
ited by an invariant Lagrangian: in this context indeed it is possible to introduce 
in a natural way and to draw a comparison between the notions of A-constant of 



- - exp(e6')w^) 
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motion and of Nocthcr A-conscrvation rule. Similarly, the symmetry properties 
of Euler-Lagrange equations and of the Hamiltonian ones can be compared, and 
some reduction techniques for the equations can be conveniently introduced. 

Finally, I have shown that the A-invariancc of the Lagrangians in the context 
of the DS is a special case of a more general and richer situation, where several 
independent variables are present and a A-conservation rule of very general form 
is true. 

Another interesting problem is the nontrivial relationship between A (or A, 
or fj,) symmetries with the standard ones. An aspect of this problem has been 
mentioned in the above section of this paper. In different situations, this may 
involve the introduction of nonlocal symmetries and other concepts in differen- 
tial geometry, as briefly indicated in the Introduction, which clearly go beyond 
the scope of the present contribution. 
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